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ABSTRACT 

With ASCA, we have detected three X-ray flares from the Class I protostar YLW15. The flares 
occurred every ~20 hours and showed an exponential decay with time constant 30-60 ks. The X-ray 
spectra are explained by a thin thermal plasma emission. The plasma temperature shows a fast-rise and 
slow-decay for each flare with kTp^ak ~ 4-6 keV. The emission measure of the plasma shows this time 
profile only for the first flare, and remains almost constant during the second and third flares at the 
level of the tail of the first flare. The peak flare luminosities Lx.peak were ~ 5-20 xlO'^^ erg s~^, which 
are among the brightest X-ray luminosities observed to date for Class I protostars. The total energy 
released in each flare was 3-6x10'^^ ergs. The first fiare is well reproduced by the quasi-static cooling 
model, which is based on solar flares, and it suggests that the plasma cools mainly radiatively, conflned 
by a semi-circular magnetic loop of length ^ 14i?0 with diameter-to-length ratio ^ 0.07. The two 
subsequent flares were consistent with the reheating of the same magnetic structure as of the first fiare. 
The large-scale magnetic structure and the periodicity of the fiares imply that the reheating events of 
the same magnetic loop originate in an interaction between the star and the disk due to the differential 
rotation. 

Subject headings: stars: flare — stars: formation — stars: individual (IRS43, YLW15) — stars: 
late-type — X-rays: spectra — X-rays: stars 



L INTRODUCTION 

Low-mass Young Stellar Objects (YSOs) evolve from 
molecular cloud cores through the protostellar (ages 
.^104-5 yr), the Classical T Tauri (CTTS: ~ 10^ yr) and 
Weak-lined T Tauri (WTTS: --10^ yr) phases to main se- 
quence. Protostars are generally associated with the Class 
and I spectral energy distributions (SEDs), which peak 
respectively in the millimeter and infrared (IR) bands. 
Bipolar flows are accompanied with this phase, suggest- 
ing dynamic gas accretion. CTTSs have still circumstellar 
disks though they have expelled or accreted the infalling 
envelopes. They are associated with the Class II spectra, 
which peak at the near-IR. Finally, as the circumstellar 
disk disappears, YSOs evolve to WTTSs, associated with 
Class III stars. Early stellar evolution is reviewed by Shu, 
Adams, & Lizano (1987) and Andre & Montmerle (1994). 

The Einstein Observatory discovered that T Tauri Stars 
(TTSs), or Class II and Class III infrared objects, are 
strong X-ray emitters, with the luminosities of 100-10000 
times larger than solar flares. These X-rays showed high 
amplitude time variability like solar flares. The temper- 
ature (~1 keV) and plasma density (rig ~10^^cm~'^) are 
comparable to those of the Sun, hence the X-ray emission 
mechanism has been thought to be a scaled-up version of 
solar X-ray emission; i.e., magnetic activity on the stellar 
surface enhanced by a dynamo process (Feigelson & De- 



Campli 1981; Montmerle et al. 1983). X-ray and other 
high energy processes in YSOs are reviewed by Feigelson 
& Montmerle (1999). 

In contrast to TTSs, Class I infrared objects are gener- 
ally surrounded by circumstellar envelopes of Ay up to ~ 
40 or more, hence are almost invisible in the optical, near 
infrared and even soft X-ray bands. The ASCA satellite, 
sensitive to high energy X-rays up to 10 keV which can 
penetrate heavy absorption, has found X-rays from Class 
I objects in the cores of the R CrA, p Oph, and Orion 
clouds at the hard band (> 2 keV) (Koyama et al. 1994; 
Koyama et al. 1996; Kamata et al. 1997, Ozawa et al. 
1999). Even in the soft X-ray band, deep exposures with 
the ROSAT Observatory detected X-rays from YLW15 in p 
Oph (Grosso et al. 1997) and CrA (Neuhauser & Preibisch 
1997). 

A notable aspect of these pioneering observations was 
the discovery of X-ray flares from Class I stars. ROSAT 
discovered a giant flare from the protostar YLW15 with to- 
tal luminosity (over the full X-ray band) of lO^'^^^^erg s~^, 
depending on the absorption. ASCA observed more details 
of X-ray flares from protostars EL29 in the Opiuchus, Rl 
in the R CrA core, and SSV63E+W in the Orion, which 
are associated with larger Nh — 10^^"^^ cm^'^ than seen 
in TTSs. 

All these findings led us to deduce that greatly enhanced 
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magnetic activity, already well-established in older TTSs, 
is present in the earlier protostellar phase. Stimulated by 
these results, and to search for further examples of the pro- 
tostellar activity in the X-ray band, we have performed 
an extensive follow-up observation of a core region in p 
Oph, with several Class I X-ray sources. The follow-up 
observation was made with ASCA 3.5 years after the first 
observation (Koyama et al. 1994, Kamata et al. 1997). 
Some previously bright Class Is became dim, while other 
Class Is were identified as hard X-ray sources. This paper 
discusses the brightest hard X-ray source, YLW15, con- 
centrating on the characteristics and implications of its 
peculiar time behavior: quasi-periodic hard X-ray flares. 
For comparison with previous results, we assume the dis- 
tance to the p Oph region to be 165 pc (Dame et al. 1987) 
throughout of this paper, although new Hipparcos data 
suggest a closer distance d ~ 120 pc (Knude & Hog 1998). 

2. OBSERVATION 

We observed the central region of p Oph cloud with 
ASCA for «100 ks on 1997 March 2-3. The telescope 
pointing coordinates were a(2000) = 16h 27.4m and 
(5(2000) = -24° 30'. All four detectors, the two Solid- 
state Imaging Spectrometers (SIS 0, SIS 1) and the two 
Gas Imaging Spectrometers (CIS 2, CIS 3) were operating 
in parallel, providing four independent data sets. Details 
of the instruments, telescope and detectors are given by 
Burke et al. (1991), Tanaka, Inoue, & Holt (1994), Ser- 
lemitsos et al. (1995), Ohashi et al. (1996), Makishima et 
al. (1996), and Gotthelf (1996). 

Each of the CIS was operated in the Pulse Height mode 
with the standard bit assignment that provides time reso- 
lutions of 62.5 ms and 0.5 s for high and medium bit rates, 
respectively. The data were post-processed to correct for 
the spatial gain non linearity. Data taken at geomagnetic 
rigidities lower than 6 GV, at elevation angles less than 
5° from the Earth, and during passage through the South 
Atlantic Anomaly were rejected. After applying these fil- 
ters, the net observing time for both GIS 2 and GIS 3 was 
94 ks. 

Each of the SIS was operated in the 4-CCD/Faint mode 
(high bit rate) and in the 2-CCD/Faint mode (medium 
bit rate). However, we concentrate on the 4-CCD/Faint 
results in this paper since YLW15 is out of the 2-CCD 
field of view. The data were corrected for spatial and gain 
non-linearity, residual dark distribution, dark frame error, 
and hot and fiickering CCD pixels using standard proce- 
dures. Data were rejected during South Atlantic Anoma- 
lies and low elevation angles as with GIS data. In order to 
avoid contamination due to light leaks through the optical 
blocking filters, we excluded data taken when the; satellite 
viewing direction was within 20° of the bright rim of the 
Earth. After applying these filters, the net observing time 
for the 4-CCD mode was 61 ks for SIS and 63 ks for SIS 
1. 

Towards the end of this observation, we detected an 
enormous fiare from T Tauri star R0Xs31 which is lo- 
cated close to YLW15 (see §3.1, source 6 in Figure 1). 
The peak flux of R0Xs31 is 1- 2 orders of magnitude larger 
than YLW15 (see §4.2), and its broad point spread func- 
tion contaminates YLW15 during the flare. Therefore, we 
excluded the GIS and SIS data taken during the flare of 
R0Xs31 in all analysis of YLW15. 



3. RESULTS AND ANALYSIS 
3.1. Images 

Figure 1 shows X-ray images of the p Ophiuchi Core F 
region in the two different energy bands (left panel: 0.7-2 
keV, right panel: 2-10 keV), obtained with SIS detectors. 
Class I sources arc indicated by crosses. Since the abso- 
lute ASCA positional errors can be as large as 40" for SISs 
(Gotthelf 1996), we compared the ASCA peak positions 
to the more accurately known IR positions of two bright 
sources in the SIS field, R0Xs21 (source 5) and R0Xs31 
(source 6), which are indicated by filled circles in Figure 
1 left panel. To obtain the ASCA peak positions, we ex- 
ecuted a two-dimensional fitting in the 0.7-2 keV band; 
we fitted these sources with a position-dependent point 
spread function in the 0.7 2 keV band and a background 
model. This procedure was done in the Display45 analysis 
software package (Ishisaki et al. 1998). The position of 
R0Xs31 was based on the flare phase of this source, while 
the position of R0Xs21 is based on the data before the 
flare of R0Xs31. IR positions are provided by Barsony et 
al. (1997). The ASCA SIS positions had an average offset 
(weighted mean by photon counts) of -1-0.18 s in right as- 
cension and —7.4" in declination from the IR frame. This 
positional offset is eorrectcxl in Figure 1. After the bore- 
sight error correction, remaining excursions between the 
X-ray and IR positions are 5.5" (rms), which is consistent 
with the SIS position uncertainty for point sources (Got- 
thelf 1996). We take the systematic positional error to be 
5.5". 

From the 2-10 keV band image, we find that the X-ray 
fluxes from Class I protostars EL29 and WL6 (sources 3 
and 4 in Figure 1, right panel) are fainter by one third 
and less than one third, respectively, comparing to those 
in the first ASCA observation made in August 1993 (Ka- 
mata et al. 1997). The brightest X-ray source in the 
2-10 keV band is an unresolved source at a(2000) = 16h 
27m 27.0s and 5(2000) = -24° 40' 50" in the position 
corrected frame. We derived this peak position by the 
two-dimensional fitting in 2-10 keV band. Since the sta- 
tistical error is 1", the overall X-ray error (including the 
systematic error) is ±6". 

The closest IR source is YLW15 with VLA position of 
a(2000) = 16h 27m 26.9s and 5(2000) = -24° 40' 49.8" 
(±0.5"; Leous et al. 1991), located 1.5" away from the 
X-ray source. The next nearest source is GY263 with IR 
position of a(2000) = 16h 27m 26.6s and 5(2000) = -24° 
40' 44.9" (±1.3"; Barsony et al. 1997). The source is lo- 
cated 5.5" from the X-ray position on the border of the X- 
ray position error circle. Thus we conclude that the hard 
X-rays are most likely due to the Class I source YLW15. 



3.2. X-ray Lightcurve of YLW15 

We extracted a lightcurve from a 3' radius circle around 
the X-ray peak of YLW15 (see Figure 1). Before the enor- 
mous flare from T Tauri star R0Xs31, which occurred in 
the last phase of this observation (see §2), we detected an- 
other large flare from Class II source SR24N, located about 
7' away from YLW15 in the GIS fleld of view. To subtract 
the time variable contamination from the SR24N in the 
extended flux of YLW15, we selected a 3' radius back- 
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0.7-2.0keV 2.0-lO.OkeV 

Fig. 1. - The SIS images of the central region of p Oph cloud in the 0.7 2 keV (left) and the 2-10 keV band (right). The contours 
are C^*O(1^0) column densities (Wilking & Lada 1983) in units of 10^*' cm^'^, from 1.0 to 2.5 in steps of 0.5, corresponding to 
Av ~30, 50, 70, 100, respectively. Class I sources are indicated by crosses (Chen et al. 1995, 1997; Motte, Andre, & Neri 1998). 
The bright sources 5 and 6 in the 0.7-2 keV band are T Tauri star ROXs21 and ROXsSl, respectively. In the 2-10 keV band 
image, the brightest source is Class I source YLW15 (source 1). Two older Class I sources, IRS51 (source 2) and EL29 (source 3), 
which have flat SEDs are also detected with a criterion of S/N > ha. EL29 and WL6 (source 4) are faint, though both of them 
emitted strong X-rays in the first ASCA observation (Koyama et al. 1994, Kamata et al. 1997). The hard X-ray source near WL6 
is a cluster of T Tauri stars (WL3, WL4, WL5; Wilking & Lada 1983). The source region and background region for YLW15 are 
also shown by circles. 



ground region (see Figure 1), equidistant from SR24N and 
YLW15. On the other hand, using such a background, we 
cannot exclude the contamination from R0Xs21 (source 5 
in Figure 1), which is 2 arcmin apart from YLW15. Since 
the X-rays from R0Xs21 are dominant below 2 keV (see 
§3.3 and Fig. 3), we examine time variability only in the 
hard X-ray band (> 2 keV) in which the flux is dominated 
by YLW15. 

Figure 2 (iippcT panel) shows the background- 
subtracted lightcurve in the 2-10 keV band with the sum 
of the SIS (SIS and 1) and GIB (GIS 2 and 3) images. 
The lightcurve shows a sawtooth pattern with three flares. 
The peak fluxes of the flares become successively less lumi- 
nous. Each flare exhibits a fast-rise and an exponential de- 
cay with an e-folding time of 31±1 ks {y^ /d.o.f. = 61/46), 
33±3 ks (80/47), and h^tlt ks (33/24), for the first, the 
second, and the third flares, respectively. We show the 
best fit lightcurves for the second and the third fiares with 
dashed lines, and show the best-fit quasi-static model (see 
§4.1.1) for the first flare with a solid line in Figure 2 upper 
panel. 
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Fig. 2. — Upper panel: Background-subtracted lightcurve of 
YLW15 obtained with all the detectors (SIS + SIS 1 + GIS 
2 + GIS 3) in the 2-10 keV band. All time bins are 1024 s 
wide. Middle panel: Time profile of the best-fit temperature. 
The error bars indicate 90% confidence. Lower panel: Same as 
the middle panel, but for the best-fit emission measure. The 
solid lines during Flare I are the best-fit quasi-static model, 
and the dashed lines during Flare II and III are the best-fit 
exponential model (see Table 2). One reheating plasma loop 
is assumed in both cases. Dotted lines during Flare II and 
III in the upper panel show each flare component without the 
contribution of the other flares, under assumption that the 
three flares occurred independently (see §4.1). Asterisks show 
the average of the model during one time bin. The time axis 
starts at 00:00:00.0 UT (1997 March 2). 
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3.3. Time-Sliced Spectra of YLW15 

To investigate the origin of the quasi-periodic flares, 
we made time-sliced spectra for the time intervals given 

in Figure 2. Wc extracted the source and background 
data for each phase from the same regions as those in the 
lightcurve analysis (see §3.2 and Figure 1). Wc foiuid that 
all the spectra show a local flux minimum at «1.2 keV. 
For example, we show the spectra obtained with SISs at 
phases 1 and 8 in Figure 3. This suggests that the spectra 
have two components, one hot and heavily absorbed, the 
other cool and less absorbed. 

Then we examined possible contamination from the 
bright, soft X-ray source R0Xs21 (source 5 in Figure 1). 
Wo extracted the spectrum of R0Xs21 from a 2' radius 
circle around its X-ray peak. We extracted the data only 
during phase 1, in order to be free from contamination 
from the flare on SR24N, which occurred during phases 
4-6 (sec §3.2). The background data for R0Xs21 were 
extracted from a 2' radius circle equidistant from YLW15 
and R0Xs21 during phase 1. After the subtraction of the 
background, the spectrum of R0Xs21 is well reproduced 
by an optically thin thermal plasma model of about 0.6 
keV temperature with absorption fixed at A'h =1.3x10^^ 
cm~^ [Ay = 0.6 mag; Bouvicr and Appenzeller 1992) with 
X^/d.o.f. = 17/23. The flux of the soft component of 
YLW15 is about 30% of the flux of R0Xs21, which is equal 
to the spill-over flux from R0Xs21. Thus the soft X-ray 
component found in YLW15 spectra is due to contamina- 
tion from the nearby bright source R0Xs21. 

Having obtained the best-fit spectrum for R0Xs21, we 
fitted the spectrum of YLW15 in each phase with a two- 
tcmpcraturc thermal plasma model. The cool component 
is set to the contamination from R0Xs21, and the hot 
component is from YLW15. For YLW15, free parameters 
are temperature (fcT), emission measure (EM), absorp- 
tion (A'h) and metal abundance. For R0Xs21, EM is the 
only free parameter and the other parameters are fixed to 
the best-fit values obtained in phase 1. We found no sig- 
nificant variation in A^h of YLW 15 from phase to phase, 
hence, we fixed the A^h to the best- fit value at phase 1. The 
resulting best-fit parameters of YLW15 for each time in- 
terval are shown in Table 1 . The best-fit spectra of phases 
1 and 8 are illustrated in Figure 3, and the time evolution 
of the best- fit parameters are shown in Figure 2. 
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The SIS (SIS -h 1) spectra of YLW15 at phase 
1 and 8 (see Figure 2). The solid lines show the best-fit two- 
temperature coronal plasma model derived from simultaneous 
fitting of GIS and SIS data. Lower panel shows the residuals 
from the best-fit model. 

4. DISCUSSION 

In this observation, we detected hard X-rays (> 2 keV) 
from the Class I source YLW15 for the first time; the 
source was not detected in the first ASCA observation ex- 
ecuted 3.5 years later (Kamata ct al. 1997). On the other 
hand, the Class I sources EL29 and WL6, which emitted 
bright hard X-rays in the first observation, became very 
faint (see §3.1). Then we conclude that hard X-rays from 
Class I protostars in the p Oph cloud are highly variable in 
the time spans, and wc suspect that the non-detection of 
hard X-rays from other Class I objects is partly due to the 
long-term time variability. Ftom YLW15, we discovered 
a peculiar time behavior: quasi-periodic flares. We shall 
now discuss the relation between each flare from YLW15 
and the physical conditions. 

4.1. Physical Parameters of the Triple Flare 

If the three intermittent flares are attributed to a single 
persistent flare, with the three events due to geometrical 
modulation, such as occupation of the flaring region by 
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Table 1 

Best-fit parameters to the spectra of YLW15°. 



Phase'' 


kT 
(keV) 


Abundance 
(solar) 


(lO^^cm-^) 


Td 

(lO^ierg s-i) 


xVd.o.f 


1 
2 
3 
4 
5 
6 
7 
8 


c c + 1.3 
^•"-1.0 

4 4+0-8 

2.9t°j 
5.3+}:? 

q + 0.4 

^••5-0.3 

1 c+0.3 
^•"-0.3 

4 2+1-^ 

q n+0.9 


o.4t°:? 
0.4^°:^ 

'-0.4 

n fi+°'i 

U-D_g 4 

U.O_Q 4 

<1.3 
<0.4 

"-"•'-0.5 


7 2+1-1 

' -^-0.8 

1+0.4 
^•J^-0.3 
9 7+0.4 
^•'-0.3 
Q 0+0.8 
"J-^-0.7 
9 q+1.5 

9 9+0.4 
^•^-0.3 

1 7+0.4 
-^•'-0.4 


15±1 
6.0±0.7 
3.4±0.8 
5.6±0.8 
4.4±1.3 
2.5±1.5 
3.5±0.5 
2.9±0.6 


181/106 
191/164 

80/84 

55/52 
163/138 

46/40 
153/139 

59/65 



Note. — The errors are at 90% confidence level. 

° Absorption column density is fixed to the best fit value for phase 1, 3.3x10^^ cm 
^ see Figure 2. 

Emission Measure, / n\dV {rie- electron density, V: emitting volume). 
Absorption-corrected X-ray luminosity in the 0.1 - 100 keV band. 



stellar rotation or orbit in an inner disk, then only the 
emission measure should have sho-wn three peaks. The 
temperature -would have smoothly decreased during the 
three fiares. Ho-wever, in our case, -we see the temperature 
foUo-wing the same pattern as the luminosity, decreasing 
after each jump, as indicated in Table 1 and Figure 2 mid- 
dle panel. To test further, -we fitted the temperatures in 
phases 1-8 with a single exponential decay model, and con- 
firmed that it is rejected -with 89 % confidence. We then 
interpret the variability of YLW15 as due to a triple fiare, 
in -which each flare foUo-wed by a cooling phase. No-w, let 
us label phases 1-3, 4-6, and 7-8 as Flare I, II, and III, 
respectively. In this subsection -we -will use the cooling 
phases to estimate the physical conditions of the plasma 
in each flare. Details are given in the Appendix. 

4.1.1. Flare I 

Here -we assume that the plasma responsible for Flare 
I is confined in one semi-circular magnetic loop -with con- 
stant cross section along the tube fiux, -with length L, and 
diameter-to-length (aspect) ratio a, based on the general 
analogy of solar-type flares. To the decay of Flare I, -we 
apply a quasi-static cooling model (van den Oord & Me-we 
1989), in -which the hot plasma cools quasi-statically as a 
result of radiative (90 %) and conductive (10 %) losses (see 
detailed comment in Appendix §1). If the cooling is truly 
quasi-static, T^^^'^ /EM sho-ws constant during the decay 
(where T is plasma temperature and EAl is the emission 
measure). We find for the three successive bins (phases 
1-3) of Flare I that {T / lO"^ K)^^ / '^f {EM/ 10^'^ cm-^) = 
61 ± 55, 59 ± 42, 25 ± 18, which are not inconsistent with 
a constant value, taking into account the large error bars. 
Then, fitting our data with the quasi-static model, we find 
satisfactory values of Xred count rates, tempera- 

ture, and emission measure (see panel 1 in Table |^). We 
conclude that our hypothesis of quasi-static cooling is well 



verified; an underlying quiescent coronal emission is not 
required to obtain a good fit. The best fits are shown by 
solid lines in Figure 2. From the peak values of T and EM, 
and the quasi-static radiative time scale, we derived loop 
parameters as listed in Table ||. The detailed procedure is 
written in Appendix. The aspect ratio a of 0.07 is within 
the range for solar active-region loops (0.06 < a© < 0.2, 
Golub et al. 1980), which supports the assumed solar anal- 
ogy- 

4.1.2. Flares II and III 

The appearance of Flare II can be interpreted either 
by the reheating of the plasma cooled during Flare I, or 
by the heating of a distinct plasma volume, independent 
from that of Flare I. In the latter case, the lightcurve and 
EM would be the sum of Flare I component and new flare 
component. The EM of the new flare component can be 
derived by subtracting the component extrapolated from 
Flare I. The derived EM are 1.6±0.4 (phase 4), 2.5±0.8 
(phase 5), and 1.9±1.5 (phase 6). Then we fitted the 
lightcurve and EM with the quasi-static model assum- 
ing the above two possibilities. However, in both of the 
cases, the model did not reproduce the lightcurve and EM 
simultaneously, and the quasi-static cooling model cannot 
be adopted throughout the triple fiare but Flare I. 

For simplicity, we fitted the parameters in Flare II with 
an exponential model. The best-fit parameters are shown 
in Table 2, and each models for the reheating and the dis- 
tinct flare assumptions are shown by the dashed and the 
dotted lines in Figure 2, respectively. The obtained Xred 
and the timescale for each lightcurve and EM were sim- 
ilar between the two assumptions. Therefore both of the 
possibilities cannot be discriminated. As for EM, both 
of the fits show no decay or very long time decay, which 
is not seen in the usual solar flares. The constant fea- 
ture in EM makes the quasi-static model unacceptable. 
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Since we cannot derive the aspect ratio of Flare II by fit- 
ting with quasi-static coohng model, assuming the ratio 
derived in Flare I and that radiative cooling is dominant, 
we deduced the plasma parameters as shown in Table 3. 
Here, we derived the values assuming the reheating sce- 
nario (the former assumption) . The results show that the 
plasma density and the volume remain roughly constant 
from Flare I. This supports that Flare II resulted from 
reheating of the plasma created by Flare I. 

As for Flare III, because of the poor statistics and short 
observed period, the fits to lightcurve with the exponential 
model give no constrain between the above two possibili- 
ties. The results are shown in Table ^ and Figure 2. We 
derived the plasma parameters of Flare III in the same way 
for Flare II, as shown in Table 3. These values are similar 
to those in Flare I and II. All these are consistent with the 
scenario that a quasi-periodic reheating makes the triple 
flare. The heating interval is '^20 hour. The loop size 
is approximately constant through the three flares and is 
as large as ~14 Rq. The periodicity and the large-scale 
magnetic structure support a scenario that an interaction 
between the star and the disk occurred by the differen- 
tial rotation and reheated the same loop periodically (e.g., 
Hayashi, Shibata, & Matsumoto 1996). The detail will be 
presented in Paper II (Montmerle et al. 1999). 

4.2. Comparison with Other Flares 

Among YSOs, TTSs have been known for strong X-ray 
time variabilities since the Einstein Observatory discov- 
ered them (see §1). At any given moment, 5-10 % are 
caught in a high-amplitude flare with timescales of hours 
(Feigelson & Montmerle 1999). The most recent example 
is that of V773 Tau, which exhibits day-long flares with 
Lx.peak ^ 2-10 xlO^^ erg s~^ and very high tempera- 
tures of ~ 10^ K (Tsuboi et al. 1998). Other examples 
of bright TTS X-ray flares are P1724, a WTTS in Orion 
{Lx,peak ~ 2 X 10^^ erg s"^; Preibisch, Neuhaiiser, & Al- 
cala'l995), and LkHa92, a CTTS in Taurus (Preibisch, 
Zinnecker, & Schmitt 1993). These X-ray properties re- 
semble those of RS CVn systems. 

Recently, a dozen protostars have been detected in X- 
rays, and four of those showed evidence for flaring (RCrA 
core; Koyama et al. 1996, EL29; Kamata et al. 1997, 
YLW15; Grosso et al. 1997, and SSV63 E-hW; Ozawa et 
al. 1999). In the "superflare" of YLW15 (Grosso et al. 
1997), an enormous X-ray luminosity was recorded during 
a few hours. If we adopt the absorption we derived {Nh — 
3x10^^ cm^^), the absorption-corrected X-ray luminosity 
in 0.1-100 keV band is Lx,peak ~ lO^'' ergs'^ The "triple 
flare" we detected in this observation does not reach the 
same level as the "superflare" : Lx,peak = 5 — 20 x lO'^^ erg 
s~^ in the same X-ray band. 

To compare the flare properties of our triple flare from 
YLW15 with other flare sources, including RS CVns, we 
selected bright flare sources as seen in Table 4. All the 
flare sources have a well-determined temperature using a 



wide range of energy band of Tenma, Ginga, and ASCA 
satellites. Since the samples of YSO flares were less, we 
added two TTS flares; the flares on R0Xs21 and SR24N 
detected in our observations (see §2 and §3). We analyzed 
them using GIS data. The densities and volumes for all 
the sources were derived assuming that radiative cooling 
is dominant. 

As a result, we found that although less energetic than 
the "superflare" , with total energies in excess of 3-6 x lO'^^ 
ergs, the triple flare are on the high end of the energy 
distribution for protostellar flares. While the plasma den- 
sities, temperatures, and then derived equipartition mag- 
netic fields are typical of stellar X-ray flares, the emitting 
volume is huge; it exceeds those in binary systems of RS 
CVns by a few orders of magnitude. 

5. SUMMARY AND CONCLUSIONS 

In the course of a long exposure of the p Oph cloud 
with ASCA, we found evidence for a 'triple flare' from the 
Class I protostar YLW15. This triple flare is the first ex- 
ample of its kind; it shows an approximate periodicity of 
~ 20 hours. Each event shows a clear decrease in the tem- 
perature, followed by reheating, with kTpeak ^ 4-6 keV, 
and luminosity Lx,peak ~ 5-20 xlO'^^ erg s~^. Apart from 
the periodicity, the characteristics of the flares are among 
the brightest X-ray detections of Class I protostars. 

A fitting with the quasi-static model (VM), which is 
based on solar flares, reproduces the flrst flare well, and 
it suggests that the plasma cools mainly radiatively, hav- 
ing semi-circular shape with length ^ 14 Rq (radius of the 
circle R ^ 4.5 Rq) and aspect ratio ~ 0.07. The minimum 
value of the confining field is ~ 150 G. The two subse- 
quent flares have weaker intensity than the first one but 
consistent with the reheating of basically the same mag- 
netic structure as the first flare. The plasma volume is 
huge; a few orders of magnitude larger than the typical 
flares in RS CVns. 

The fact that the X-ray flaring is periodic suggests that 
the cause for the heating is periodic, hence is linked with 
rotation in the inner parts of the protostar. The large size 
of the magnetic structure and the periodicity support the 
scenario that the flaring episode has originated in a star- 
disk interaction; differential rotation between the star and 
disk would amplify and release the magnetic energy in one 
rotation period or less, reheating the flare loop as observed 
in the second and third flares. 
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Table 2 

Fitting of cooling phase of flares 



TTlOrlpl P'VTiFP'^dlOTl 


T 


ol 




'J- 


y2 
Area 














[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


^I,qsW = X {1 + - tl)/3T}-« 


J [cntsks"-'-] 


4 


172 ±5 


31 ±1 


0.8 (45) 




kT [kcV] 


8/7 


5.6 ±0.6 


[31] 


0.5 (2) 




EM \\{f^ cni-3] 


26/7 


10.4 ±0.9 


[31] 


0.4 (2) 


J'llit) =:FpX exp{-{t - tu)/T} 


/ [ciitsks^^] 




62 ±3 


33 ±3 


1.7 (47) 




fcr [kcV] 




4± 1 


50 ±20 


1.5 (1) 




EM {Wf'^ cni-3] 




3± 1 


> 53 


0.4 (1) 




/ [cntsks^i] 




46 ± 3 


29 ±2 


1.6 (47) 




EM [10" ^;m--'] 




l.G ± 0.8 


> 54 


U.(i (1) 


^iii(i) = X exp{-{t - iiii)/r} 


/ [ciitsks^-'-] 




35 ±2 


60 ±20 


1.4 (24) 




kT [kcV] 




5±3 


55 (> 16) 


= (0) 




EM [10-'^4 cm--^] 




2.5 ±0.7 


71 (> 23) 


= (0) 




I [cntsks"^] 




28 ±2 


60 ±20 


1.4 (24) 



Note. — Panel 1 gives the quasi-static fit for Flare I. Panel 2 (Panel 3) gives the exponential fit for Flare II (Flare III): 
first assuming that the flares are reheating of the same plasma volume as the previous flare; second assuming that the 
three flares are independent (the previous flares contribute to the latter flares). t\, tn, and t\\\ are the beginning time of 
Flare I, II, and III, respectively. We assume only positive values for the decay time (Column [5]). The errors are at 90 % 
confidence level. 



Table 3 

Physical parameters for the triple flare 



Flare 


a 


L 


rie 


B 


Lx.peak 




Etot 






[Re] 


[10" cm-3] 


[G] 


[10^1 erg s-1] 


[ks] 


[10^6 erg] 


I 


0.07±0.02 


14±2 


0.5±0.1 


150±20 


20±1 


31±1 


6.0±0.5 


II 


[0.07±0.02] 


11±5 


0.4±0.1 


120±30 


8±1 


33±3 


2.8±0.6 


III 


[0.07±0.02] 


15±10 


0.3±0.2 


100±60 


4.6±0.5 


60±20 


2.7±1.0 



Note. — a: diamctcr-to-lcngth ratio of the loop; L: loop length; rig: electron density; B: minimum value of the mag- 
netic field confining the loop plasma; Lx,peak'- X-ray luminosity at the flare peak in the 0.1-100 keV band; t: characteristic 
decay time of the lightcurve; Etot'- total energy released during the cooling phase in X-ray wavelength. 
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YLW15 



Table 4 

Comparison with the other flares 



Source 


SED 




V 


T 




EMpeak 


Etot 


ref. 






[10" cin-3] 


[10^3 cm^] 


[ks] 


[keV] 


[10^^ cm-3] 


[10^6 erg] 




YLW15(I) 


I 


0.5 


3 


31 


6 


0.7 


6 




YLW15(II) 


I 


0.4 


2 


33 


5 


0.3 


3 




YLW15(III) 


I 


0.3 


2 


60 


4 


0.2 


3 




EL29 


I 


1 


0.2 


8 


4 


0.3 


0.4 


(1) 


RCrA 


I 


1 


0.04 


20 


6 


0.04 


0.2 


(2) 


SSV63E+W 


I 


1 


2 


10 


6 


2 


4 


(3) 


SR24N 


II 


5 


0.1 


6 


7 


3 


6 




ROXsSl 


III 


3 


0.6 


7 


6 


6 


13 




R0X7 


III 


3 


0.01 


6 


3 


0.1 


0.1 


(1) 


V773Tau 


III 


3 


0.6 


8 


10 


5 


8 


(4) 


Algol 




1 


0.09 


20 


5 


0.1 


0.6 


(5) 


TT Peg 




8 


0.006 


2 


6 


0.3 


0.2 


(6) 


Au :\iic 




() 


0.002 


2 


5 


O.OG 


o.o.s 


(7) 



Note. reference; (1) Kamata et al. 1997 

(2) Koyama ct al. 1996 and Hamaguchi private comm. 

(3) Ozawa et al. 1999 

(4) Tsuboi et al. 1998 

(5) Stern ct al. 1992 

(6) Doyle et al. 1991 

(7) Cully et al. 1993 
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APPENDIX 
APPENDIX 

DETERMINATION OF THE PHYSICAL PARAMETERS OF THE FLARES 
1. LOOP PARAMETERS 

We make use of the general treatment for solar-type flares, put forward by van den Oord & Mewe (1989, hereafter 
VM). The decrease in the thermal energy of the cooling plasma is assumed to be caused by radiative (rr) and conductive 
losses (tc): its decay time is thus 1/tcs = I/ti- + 1/tc. VM assumed that the flare Hghtcurve and temperature decrease 
exponentially with decay times and tt, respectively. This effective cooling time is related to observed time scales by 
^/tcS = 7/8tt + l/2Td. We assume here that the flare occurs in only one semi-circular loop with constant section along 
the tube flux, radius R {R = L/tt; length L), diameter-to-length ratio a, and volume V. VM gives an expression of L 
versus a, depending on TcS, the temperature, the emission measure (hereafter EM), and the ratio Tr/rc. Because of the 
assumed exponential behavior of the lightcurve and temperature, the moment at which this expression is applied is not 
important. The only restriction is that both the temperature and the EM started to decrease. 

Due to the low statistics, we have only time-sliced values of the temperature and the EM. Let us call ti (tf) the 
beginning (end) of the time interval within which the temperature was estimated, we have: T = /;/ T{t') dt'/{tf - U). 
We use this relation to find the behavior of the temperature as a function of time. 

We now have a relation between L and a, but the ratio t-c/t,, is unknown, and even worse may change during the 
decay of the flare. An exception to this is when the flare volume cools quasi-statically, evolving through a sequence of 
quasi-static equilibria, where the loop satisfies scaling laws, and where T]./tc — est. Due to the dependency of Tj. and Tc 
on the temperature and the EM, t^/tc — est implies T^^/^ /EM — est. VM gives in that case analytical expressions for 
several physical quantities versus time, all depending on the quasi-static radiative time scale (Tr,qs). Tr,qs can be estimated 
from the lightcurve which must be proportional to the radiative loss. The temperature and the EM can be fitted as 
described above using this value of rr,qs, and give the peak values kT^^qs, EMp^qs (see Col. [4]-[5] in our Table ||, panel 1 
for details). 

Using an expression for the radiative time (eq.[23a] of VM), the EM (eq.[3] of VM) and the scaling law (SL[2], eq.[196] 



of VM), we obtained the loop characteristics for the quasi-static model: 

a = 1.38 X (r,,qs/10 ks)-^^^ x {kTp^qJkeV)-^^/^^ x (i^A/p.qs/lO^" em~^)^/^, (1) 

L - 1.0 i?0 X (r,,qs/10 ks) X {kTp^q,/keVy^^ (2) 

He = 4.4 X 10^" em-^ x (r^.q^/lO ksy^ x {kTp^q,/keVf/^ . (3) 



Using an expression for the ratio t-^/tc on page 252 of VM, we found Ti /tc = ji x f, with the parameter fi depending 
only on the exponent of the temperature in the expression for the radiative loss (for temperature above 20 MK, fj, = 0.18), 
and the multiplicative function / coming from the expression for the mean conductive energy loss (formula [7] of VM) , 
which is equal to 4/7 for a loop with a constant section. Thus, t^/tc = 0.1. ^ In other words, the quasi-static model 
implies that 91 % of the lost energy are radiative losses: radiation is the dominant energy loss process. 

Assuming that the cooling is only radiative we used these simplified relations based on the exponential decay of the 
lightcurve, the temperature, and the EM: 

Tie = 4.4 X 10^" X (Td/lO/cs)-^ X (/^Tp/fceT/)3/^ for kT > 2keV, (4) 

L = 7.4 i?0 X (a/0.07)-2/3 x (rd/lO /cs)^/^ x [kT^/keVy^'^ x (£;A^p/10^'')l/^ (5) 
with fcTp (EMp) the peak value of the temperature (EM). 

2. MAGNETIC FIELD 

Assuming equipartition between the magnetic pressure B'^/Stt and the ionized gas pressure 2n(.kT, we can obtain a 
minimum value of the magnetic field confining the emitting plasma using: 

B = 28.4 G X (ne/10^" em-^f'^ x [kT/keV f^'^ (6) 
3. RELEASED ENERGY 

For estimating the energy released by the flare during its cooling phase we need the peak luminosity value of this flare. 
As the lightcurve must be proportional to the intrinsic luminosity, we flt the time averaged intrinsic luminosities in the 
0.1-100 keV band (given in Table 1) with the same model used for the lightcurve fitting. This gives the peak luminosity, 
-f'X.peak, and a characteristic decay time r. Thus, the total energy released by this flare is: 

E,,t - 10^' erg x (Lx, peak/ W erg s'^) x (r/fcs) (7) 



^VM wrote Tr/rc = fi = 0.18, and the analytical expression for the conductive energy loss in the quasi-static model without taking this 
multiplicative factor 4/7 into account (see Table 5 of VM). 



